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Abstract. The idea of dissipative mechanical system with delay is 
proposed. The paper studies the phenomenon of dissipation with delay 
for Euler-Poincare systems on Lie algebras or equivalently, for Lie- 
Poisson systems on the duals of Lie algebras. The study was suggested 
by the work [2] and it is ended with a discussion regarding the stability 
and the Hopf bifurcations for the free rigid body with delay. 
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Introduction 

In many applications one assumes the system under consideration is 
governed by a principle of causality; that is, the future state of the system 
is independent of the past states and is determined solely by the present. 
If it is also assumed that the system is governed by an equation involving 
the state and the rate of change of the state, then, generally, one considers 
either ordinary or partial differential equations. However, under a closer 
scrutiny, it becomes apparent that the principle of causality is often only 
a first approximation to the true situation and that a more realistic model 
would include some of the past states of the system. The simplest type 
of past dependence is through the state variable and not the derivative of 
the state variable, the so-called retarded functional differential equations or 
retarded difference equations or systems with delay. Systems with delay are 
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studied in many biological research topics, as well as in several branches of 
engineering, in neural networks, in economics, optimal production decision 
for an oligopoly with information lags ([9], [11], [12] etc.). 

Functional differential equations on finite-dimensional manifolds are con- 
sidered in [5], where the topological properties of the global attractor of an 
retarded functional differential equation in terms of limit capacity and Haus- 
dorff dimension are presented. 

The purpose of this paper is to study the phenomenon of dissipation 
with delay for Euler-Poincare systems on Lie algebras or equivalently, for 
Lie-Poisson systems on the duals of Lie algebras. The dissipation without 
delay inducing instabilities for Euler-Poincare systems is studied in [2]. The 
dissipation with delay that we construct has the essential feature: the energy 
is dissipated, but the angular momentum is not. In the context of Euler- 
Poincare or Lic-Poisson systems this means that the coadjoint orbits remain 
invariant but on them the energy is decreasing along orbits. It is interesting 
the geometry behind the construction of the nonlinear dissipative term with 
delay, which has a Brockett with delay double bracket form. In fact, this 
form is well adapted to the study of dissipation with delay on Lie groups 
since it was constructed as a gradient system and it is well known in other 
contexts that this formalism plays an important role in the study of integrable 
systems. 

The general equations of motion for dissipative systems with delay that 
we consider have the following form: 



where H is the total energy of the system, {F, H} is a skew symmetric 
bracket which is a Poisson bracket in the usual sense and where | |F, if 1 1 



is a symmetric bracket with delay. 

The outline of this paper is as follows. In section 1 the functional diffe- 
rential equations on manifolds are introduced. In section 2 a concrete idea of 
dissipative mechanism with delay is presented. The same formalism can be 
applied to other systems as well. Some of the basic and essential facts about 
dissipative mechanical systems with delay are described in section 3. In sec- 
tion 4 we study Lagrangian systems that are invariant under a group action 
and we shall add to them, in the sense of section 3, dissipative fields with 
delay that are equi variant. It is given by a necessary and sufficient condition 



f = {f,h}-{{f,h}}^ 




2 



that the integral curves of the vector field Z + Y, for a vertical G-invariant 
vector field Y on TQ x TQ and a Lagrangian vector field Z of a G-invariant 
Lagrange function, preserve the inverse images of the coadjoint orbits in g* 
by the momentum map J. In section 5 the Eulcr-Poincare and Lie-Poisson 
equations with delay are given and the double bracket with delay is defined; 
in the case "without delay" this is defined in [2]. In section 6 is presented 
the free rigid body with delay if the components (Ji, I2, 13) of the moment of 
the inertial tensor satisfy the conditions Ii > I2, h > h- For the equilibrium 
state ill = (m//i,0,0)^, m ^ 0, the value Tc of the delay for which fli is 
asymptoticaly stable and the value tq of the delay for which there is a Hopf 
bifurcation arc determined. They are also determined the local center man- 
ifold and the quantities Gi(0), /i2,72,/32 sketching out the direction of the 
Hopf bifurcation, the stability and the period of the bifurcating periodic so- 
lutions. These quantities are calculated for fixed values of a and m. Finally, 
some conclusions are drawn and further research directions are discussed in 
the last section. 



1. Functional differential equations on manifolds 

In this section, we begin with examples that will serve as a motivation 
for the consideration of functional differential equations on manifolds. 
Example 1. For any constant c, the scalar equation 

q{t)^csm{q{t-l)) (1.1) 

can be considered as an retarded functional differential equations (RFDE) 
on the circle S ^ = [x e IR^ {x^y + (x^)^ = l| by considering q as an angle 
variable only determined up to a multiple of 27r. 

Example 2. If b, c are constants then we write the second-order RFDE 

q{t) + bq{t) ^csm{q{t-l)) (1.2) 

as a system of first-order RFDE 

eit)=q\t), 

q^{t) = csm{q\t - 1)) - bq^{t). ^ ' ' 

By considering as an angle variable only determined up to a multiple 27r, 
the equation (1.3) is a RFDE on the cylinder x IR. We remark that 
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we can take the space of initial data for the solution {q^,q^) of (1.3) as 
C([-l,0],5^) X H. 

The simplest type of past dependence in a differential equation is that in 
which the past dependence is through the state variable and not the derivative 
of the state variable, the so - called retarded functional differential equations 
or retarded differential difference equations. For a discussion of the physical 
applications of the differential difference equation 

q{t)^F{t,q{t),q{t-T), 

to control problems, see [5], [6]. 

Example 3. The equations involved in the study of vibrating masses 
attached to an elastic bar are 

qHt)+co!qHt)=eMq\t),q' (t) ,q'{t),e (t) )+liq\t-T), 
e{t)+ulq\t)^eh{q\t)A\t),q\t)A\t))+^2q\t-T). ^ ' ' 

Example 4. Let S'' = {{q\ q^ q^) e 1R^ {q^f + {q^f + {q'f = l} and 
consider the following system of RFDE: 

q\t) = ~q\t-r)q\t)-q\t), 

q\t) = q\t-T)q\t)-q\t), (1.5) 
q\t) = q\t) + q\t). 

If {q^{t), q^{t), q^{t)) is a solution of the equation (1.5), it is easy to see that 

j:q\tmt) = ^ 

1=1 

for all i > 0, Vr. As a consequence, for t > 0, (9*(^)) = , a constant. 

i=l 

Thus, if an initial condition = {(f)^,(f)'^,(f)^) satisfies (j){6) G S"^ for all 9 G 
[— r, 0], we conclude that the solution {q^,q^,q^)(t; (/)) G for all t > 0. 
With this remark, we can define an RFDE on S"^ by the map 

F : (/) = (0\ </)^ </)3) G C{[-T, 0], S') ^ F(0) G TS' 

where F((f)) is the tangent vector to S'^ at the point 0(0) given by 

F(0) = (-0i(-r)02(O) - 0^(0), <f>\-T)<f>\0) - 0=^(0), 0^(0) + 0^(0)). 
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We now formalize the notions in these examples to obtain a RFDE on a 
n-dimensional manifold. Roughly speaking, a RFDE on a manifold Q is a 
function F mapping each continuous path lying on Q, e C([— r, 0], Q], 
into a vector tangent F{(f)) to Q at the point 0(0) e Q. 

Let Q be a separable finite n-dimensional manifold (configuration 
manifold), / = [— r, 0], r > 0, and C{I,Q) the totality of the continuous 
maps (j) of I into Q. The space C{I,Q) is separable and is a C°°-manifold 
modeled on a separable Banach space. 

If p : C{I,Q) Q is the evaluation map, p(0) = 0(0), then p is C°°, 
and for each q & Q, p~^{q) is a closed submanifold of C{I, Q) of codimension 
n = dimQ. 

A retarded functional differential equation (RFDE) on Q is a continuous 
function F : C{I ,Q) ^ TQ such that tttq o F — p. If we want to emphasize 
the function F defining the RFDE, we write RFDE(F). 

A solution of RFDE(F) is defined in the obvious way, namely, as a conti- 
nuous function q : [— r, a) Q, a > 0, such that q{t) exists and is continuous 
for t e [0,a) and {q{t),q{t)) = F{qt), for t G [0, a) where qt{9) = q{t + 6), 
9 e [— T, 0]. Locally, if -F(0) = (0(0), /(0)) for an appropiate function /, 
then this is equivalent to q{t) — f{qt)- 

The basis theory of existence, uniqueness and continuous dependence on 
initial data for general RFDE on manifold is the same as the theory when 

g = iR" [4]. 

Any C'^-vector field on Q defines a C'^-RFDE on Q. In fact, \l X : Q ^ 
TQ is a C'^-vector field on Q, it is easy to see that F = X o p is a C'^-RFDE 
on Q. 

To show that the equation considered in Example 2 is an RFDE according 
to the definition, we need the concept of a second order RFDE on Q. Let 
F : C{I,TQ) — > TQ x TQ be a continuous function that locally has the 
representation 

F(0, V) = ((0(0), v(o)) , (v(o), /(0, m ■ 

The solution {x{t),y{t)) of the RFDE(F) on TQ satisfies the equations 

Q{t)=y{t), y{t) = f{xt,yt) 

where q{t) G Q. If it is possible to perform the differentiations, then we 
obtain the second-order equation 

m = fi<lt,qt)- (1.6) 
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If we now return to Example 2, we sec that the formulation requires that 
we consider initial data in the space C(/, S^) x C{I, IR). However, this does 
not affect the dynamics since the solution will be in space C(/, S^) x H after 
one unit of time. 

Let g : QxQ ^ TQ be such that {7iTQ°g){(l, q) = Q and let d : C{I,Q) ^ 



QxQhe defined by = (0(0), 0(-r)). The function F = god is a RFDE 
on Q which can be written locally as 



where g{q, q) = {q, g{q, q)). The equations (1.7) is a delay differential equation 
(DDE) on Q. 

A C^-vector field with delay is given by the mappings g : Q x Q ^ T*Q, 
such that {TTT*Qg){qi, ^2) = 9i, and d : C(/, Q) ^ Q x Q, defined by d{(j)) = 
(0(0), 0(— r)). The function G = gr o d is a covector field with delay on Q 
that can be written locally as 



where g{q,q) = {q,f{q,q)). 

The topological properties of the RFDE on manifolds are discussed in [5] . 

2. Motivating examples 

To get a concrete idea of the type of dissipative mechanisms with delay wc 
have in mind, we now give a simple example of it for perhaps the most basic of 
Euler-Poincare or Lie-Poisson systems, namely the rigid body. Here, the Lie 
algebra in question is that of the rotation group; that is, the Euclidean space 
IR^ interpreted as the space of body angular velocities Q equipped with the 
cross product as the Lie bracket. On this space, we put the standard kinetic 

energy Lagrangian L(Q) = -(/ • Q) [where / = diag (/i, /2, -^3) is the moment 

of incrtial tensor] so that the general Euler-Poincare equations become the 
standard rigid body equations for a freely spinning rigid body: 



q{t) = g{q^),q{t-T)) 



(1-7) 



u^{t)^h{q{t),q{t-T))dq' 



(1-8) 



m = (7Q) X Q 
or, in terms of the body angular momentum M — IQ,, 



(2.1) 



M = MxVl. 



(2.2) 
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In this case, the energy equals the Lagrangian; E{Q) = L{Q) and the energy 
is conserved by the solutions of (2.1). Now we modify the equations by 
adding a term with delay 

M = M xfl + aM X (^M xn^ (2.3) 

where a is a positive constant, M{t) — M{t — r), for i > 0, r > and 
M{t) = (fit), for t e [-T,0]. 

A related example is the Laudau-Lifschitz equations with delay for the 
magnetization vector M in a given magnetic field B 

M ^-fM X B + - — (m X (m X b)) (2.4) 

' ||Mp-cos^^ \ JJ ^ ^ 

where 7 is the magneto-mechanical ratio, A is the damping coefficient due 
to domain walls and 9 is the angle between M and M. 

One checks in each case that the addition of the dissipative term with 
delay has a member of interesting properties. First of all, this dissipation 
with delay is derivable from a 5*0 (3)- invariant force field. However, it is 
induced by a dissipation function with delay in the following restricted sense. 
It is a gradient when is restricted to each momentum sphere (coadjoint orbit) 
where each sphere carries a special metric (later to be called the normal 
metric). Namely, the extra dissipative term with delay in (2.3) equals the 
negative gradient of the Hamiltonian with respect to the following metric 
on the sphere. Each vector v e JR^ can be orthogonally decomposed with 
respect to the standard metric on IR^ into^a component tangent at M to the 
sphere ||M|P = and a component on M, where M 

M ■ V 



-M- 



cos 6 cos 6 



M x(m X v)\ (2.5) 



where 9 is the angle from M and M. The metric on the sphere is chose to be 
(c^cos^)~'^q; times the standard inner product of the components tangent to 
the sphere in the case of the rigid body model with delay and just A times 
the standard metric in the case of the Laudau-Lifschitz equations with delay. 

Secondly, the dissipation with delay to the equations has the obvious form 
of a repeated Lie bracket, i.e. a double bracket, and it has the properties 
that the conservation law 

|||Mf = (2.6) 
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is preserved by the dissipation with delay (since the extra force is orthogonal 
to M) and the energy is strictly monotone except at relative equilibria. In 
fact, we have: 

^ " M X Q ^ (2.7) 



dt 



for the rigid body and 



d 
'dt 



E 



-a 



A 



||M||2cos^ 



M X 5 



(2.8) 



in the case of the Laudau-Lifschitz equations, so that the trajectories on the 
angular momentum sphere converge to the minimum (for a and A positive) 
of the energy restricted to the sphere, apart from the set of measure zero 
consisting of orbits that are relative equilibria or are the stable manifolds of 
the perturbed saddle point. 

Another interesting feature of the dissipations with delay is that they can 
be derived from a bracket in the same way that the Hamiltonian equations 
can be derived from a skew symmetric Poisson bracket. For the case of the 
rigid body with delay, this bracket is 

{{F, K}} = a{M X VF) • (m x VK) . 

As we have already indicated, the same formalism can be applied to other 
systems as well. In fact, later in the paper we develop an abstract construc- 
tions for dissipative with delay terms with the same general properties as the 
above examples. 



3. Dissipative systems with delay 

For later use, it will be useful some of the basic and essential facts about 
dissipative mechanical systems with delay. Let Q be a manifold, L : TQ IR 
be a smooth function and let tt : TQ ^ Q he the tangent bundle projection. 
Let WL : TQ — > T*Q be the fibre derivative of L; recall that it is defined by 



< WL{v),w > 



d_ 

de 



e=0 



L{v + ew). 



(3.1) 



Where <, > denotes the pairing between the tangent and cotangent spaces. 
We also recall that the vertical lift of a vector w G TgQ along v e TqQ is 
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defined by: 




{v + ew) e Ty{TQ). 



(3.2) 



The action and energy of L are defined by 



A{v) 



< FL{v),v > 



(3.3) 



and 



A{v)~L{v). 



(3.4) 



Let = (]FL)*0 denote the pull back of the canonical sympletic form 
on T*Q by the fibre derivative of L. 

A vector field Z on TQ is called a Lagrangian vector field of L if 



In this generality, Z need not exist, nor be unique. However, we shall assume 
throughout that Z is a second order cqiiation; that is Tn o Z is the identity 
on TQ. A second order equation is a Lagrangian vector field if and only if 
the Euler-Lagrange equations hold in local charts. We note that, by skew 
symmetry of Ql, the energy is always conserved; that is, E is constant along 
an integral curve of Z. We also recall that the Lagrangian is called regular 
if Ql is a (weak) sympletic form; it is nondcgcncratc. This is equivalent 
to the second fibre derivative of the Lagrangian being, in local charts, also 
weakly nondcgcncratc. In the regular case, if the Lagrangian vector field 
exists, it is unique, and is given by the Hamiltonian vector field with energy 
E relative to the sympletic form Q^. If, in addition, the fibre derivative 
is a global diffeomorphism, then Z is the pull back by the fibre derivative 
of the Hamiltonian vector field on the cotangent bundle with Hamiltonian 
H = E o (JFL)^^. It is well known how one can pass back and forth between 
the Hamiltonian and Lagrangian pictures in the hyperregular case [2]. 

Consider a general Lagrangian vector field Z for a (not necessarily regular) 
Lagrangian on TQ. A map Y : TQ x TQ — > TQ is called a dissipative vector 
field with delay if is vertical, i.e. TnoY — and if at each point of TQ x TQ 



If the inequality is pointwise strict at each nonzero v G TQ, v G TQ, v ^ v, 
then we say that the map Y is dissipative. A dissipative Lagrangian system 



iz^L — dE. 



(3.5) 



<dE,Y> < 0. 



(3.6) 
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on TQ is a vector field of tlie form X = Z + Y , where Z is a (second order) 
Lagrangian vector field and F is a dissipative vector field with delay. Define 
the 1-form on TQ x TQ by 

= -iy^L (3.7) 
and the force field with delay -.TQxTQ^ T*Q given by 

< {v, v),w>= A'^iv, v)-V, = -n^iv) {Y {v, v) , K) (3.8) 

where TniV-u) = w, and K. G T>v{TQ). 

Proposition 3.1. A vertical vector field Y : TQ x TQ TQ is dissi- 
pative with delay if and only if the induced force field with delay F^ satisfies 
< F^ {v,v) ,v > < for all nonzero v e TQ (< for the weakly dissipative 
with delay). 

Proof. Let F be a vertical vector field Y : TQ x TQ ^ TQ with 
Ttt oY = 0, A^ the form on TQ x TQ given by (3.7) and the force field 
with delay given by (3.8). If Z denotes the Lagrangian system defined by L, 
we get 

{dE-Y){v,v) = {tznL){Y){v,v) = QL{Z,Y){v,v)^ 
= -nL{v){Y{v,v),Z{v)) = 
= < F"" {v,v),T,7r{Z{v)) >=< F^ {v,v) ,v >, 

since Z is a second-order equation. We conclude that < dE, F > < if and 
only if < {v,v) ,v > < 0, for all {v, v) e TQ x TQ, v ^ v, which gives the 
result. 

Treating A^ as the exterior force with delay 1-form acting on a mechani- 
cal system with a Lagrangian L, we now shall write the governing equation of 
motion. The basic principle is of course the Lagrange-d'Alembert principle. 

The Lagrangian force associated with a given Lagrangian L and a given 
second-order vector field X is the horizontal 1-form on TQ defined by 

MX) = ix^L - dE. (3.9) 

Given a horizontal 1-form cu (referred to 1-form as the exterior force with 
delay), the local Lagrange-d'Alembert principle states that 

0L(X)+a; = O. (3.10) 
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It is easy to check that is indeed horizontal if X is of second order. 

Conversely, if L is regular and if is horizontal, then X is of second 

order. One can also formulate an equivalent principle in variational form. 

Given a Lagrangian L and a force field with delay (as defined in Proposi- 
tion 1) the integral Lagrange-d'Alembert principle with delay for a curve q{t) 
in Q is 

S f L {q{t), m) dt+ ['f ({q{t), q{t)) , ^(*))) • Sqdt = 0, (3.11) 

where the variation is given by the usual expression 

^/;Lfaw,,w)*^/;(|-|(|))^,-* (3.12) 

for a given variation 5q (vanishing at the endpoints). 

The two forms of the Lagrange-d'Alembert principle are equivalent. This 
follows from the fact that both give the Euler-Lagrange equations with for- 
cing with delay in local coordinates (provided that Z is of second order). We 
shall see this in the following development. 

Proposition 3.2. Let ou be the delay exterior force 1-form associated to 
a vertical vector field with delay Y , i.e. lo = = —iy^l- Then X = Z + Y 
satisfies the local Lagrange-d'Alembert principle with delay. Conversely, if, 
in addition, L is regular, the only second-order vector field X satisfying the 
local Lagrange-d'Alembert principle with delay is X — Z -\-Y. 

Proof. For the first part, the equahty 0l(^) -|- a; = is a simple ve- 
rification. For the converse, we already know that X is a solution and the 
uniqueness is guaranteed by regularity. 

To develop the differential equations associated to X = Z -\- Y, we 
take uj = —iy^L and note that, in a coordinate chart, Y ({q,q) (qiQ)^ = 
d 

((?) q) ) ?)) 'K~ and the equation (3.10) is given by 



dq^ 

d (dL\ dL d'^L 



r - 



dt \dq^J dq^ dcf-dq- 
The force 1-form with delay is therefore given by 

d^L 



Y^({q\q'),(q\it)). (3.13) 



{q\q% {q\'^)) = ^^Y^ {{q\q'), {q',^')) dq' (3.14) 
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and the corresponding force field with delay is given by 

Thus, the condition for an integral curve takes the form of the Euler- 
Lagrange equations with force with delay 

Since the integral Lagrange-d'Alembert principle with delay gives the 
same equations, it follows that the two principles are equivalent. From now 
on, we shall refer to either one as simple the Lagrange-d'Almbert principle 
with delay. 

Example 1. The inertial neuron with time delay is described by the 
equations 

n 

q^ = -aq^-bq^ + cf{q'-h^) + dJ2fiQ'-hq'), ^ = (3.17) 
where a, 6, c, d > 0, /i > are constants and q^{t) — q^{t — T) is the time delay. 

1 " 2 " 

For n — 1 the model is discussed in [1]. For L(q,q) — (^g*j — 

2 i=i i=i 

and the force field with delay given by 

Y' {q, q) = -aq' + cf(q' - h^) + dj^ f (q' - hq^) , (3-18) 

the equations (3.17) are the Euler-Lagrange equations with force with delay. 
The model for n = 2 is analyzed in [9]. The force field with delay (3.18) is 
not dissipative. 

Example 2. The simplest mechanical model of the regenerative machine 
tool vibration in the case of the so-called orthogonal cutting is given by the 
equation 

q + 2kaq + a^q = —f{q,l3), (3.19) 
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where a = ys/m is the natural angular frequency of the undamped free 
oscillating system, k is the so-called relative damping factor and / {q, (3) is 
the cutting force as a function of technological parameters and as a function 
of the chip thickness /3 which depends on the position q of the tool edge. In 
[12] the function / is given by 



qiY 



(3.20) 



where ki — -awP 4 is the parameter depending on further technological 

parameters and considered to be constant, w is the width of the chip, Q 
is the constant angular velocity of rotating work-piece (or tool), the delay 

T = — is the time period of one revolution and q^{t) — q^{t — t). For 
L (?) q) — 7^0^ ~ o^;?^ the force filed with delay given by 



Y{q, q)^-2kaq + f{q, q, (5) , 



(3.21) 



the equation (3.19) is the Eulcr-Lagrange equation with force with delay. The 
model was analysed in [12]. The force field given by (3.21) is not dissipative. 



4. Equivariant dissipation with delay 

In this section we study Lagrangian systems that are invariant under a 
group action and we shall add to them, in the sense of the preceding sections, 
dissipative fields with delay that are equivariant. This invariance property 
will yield dissipative mechanisms with delay that preserve the basic conserved 
quantities, yet dissipate enery, as we shall see. 

Let G be a Lie group acting on the configuration manifold Q and assume 
that the lifted action leaves the Lagrangian L invariant. In this case, the 
fibre derivative WL : TQ T*Q is equivariant with respect to this action 
on TQ and the dual action on T*Q. The action A, the energy E and the 
Lagrangian 2-form Q,l are all invariant under the action of G on TQ. Let Z 
be the Lagrangian vector field for the Lagrangian L, which we assume to be 
regular. Because of regularity, the vector field Z is also invariant under G. If 
the action is free and proper, so that (TQ) /G is a manifold, then the vector 
field and its fiow Ft drop to a vector field Z'^ and a fiow on {TQ)/G, [2]. 



13 



Let J : TQ — > g* be the momentum map associated with the G-action, 
given by: 

J{v,)-^^<WL{v,),^Q{q)> (4.1) 

for Vq e TqQ and for ^ G g, where denotes the infinitesimal generator for 
the action on Q. The infinitesimal generator for the action on the tangent 
bundle will be likewise denoted by ^tq and for later use, we note the relation 
TttoC^tq = ^q°it- If v{t) denotes an integral curve of the vector field with an 
equivariant dissipation term with delay Y added as in the preceding section 
and if J^{v) —< J{v),^ > is the ^-component of the momentum mapping, 
then we have: 

j/{v{t)) = dJHv{t)) ■ Z{v{t)) + dJHv{t))Y {v{t),v{t)) ■ (4.2) 

The first term vanishes by conservation of the momentum map for the La- 
grangian vector field Z. Prom (2.8) and the definition of the momentum 
map, we get: 

dJ^v) ■ Y {v,v) = {ti,^nL){Y) {v,v) = -{tyn^){^rQ) {v,v) 

= <F''{v,v),T,T7:{^tq(v))>^ (4.3) 
= <F''iv,v),CQ{7r{v))> 

and therefore 

pHv{t)) =< F''{v{t),v{t)),^QTn{v{t)) > . (4.4) 

We summarize this discussion as follows. 

Proposition 4.1. The momentum map J : TQ — >• g* is conserved under 
the flow of a G -invariant dissipative vector field with delay Z + Y if and only 
if < , o T >— for all Lie algebra elements ^ G g. 

In this paper we shall consider dissipative vector fields with delay for 
which the flow drops to the reduced spaces. Thus a first requirement is that 
y be a vertical G-invariant vector field on TQ x TQ. A second requirement 
is that all integral curves v{t) oi Z + Y preserve the sets J~^& , where & is 
an arbitrary coadjoint orbit in g*. Under these hypotheses the vector field 
Z + y induces a vector field Z^ + Y^ on {TQ x TQ)/G x G that preserves 
the symplectic leaves of this Poisson manifold, namely all reduced spaces 
J-^&/G. 
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The condition that v{t) G J ^O" is equivalent to J{v(t)) G i.e., to 
the existence of an element ri{t) = ri{v{t)) G g such that dJ{v{t)) / dt = 

dJlMR = jlvm](^y(^t)) (4.5) 
dt 

for all C ^ S- III view of (4.4), we get the following 

Corollary 4.2. The integral curves of the vector field Z + Y , for Y a 
vertical G -invariant vector field on TQ x TQ and Z the Lagrangian vector 
field of a G-invariant Lagrangian function L : TQ — > TR, preserve the inverse 
images of the coadjoint orbits in g* by the momentum map J if and only if 
for each {v,v) G TQ x TQ there is some r){v,v) G g such that 

< {v,v) , {^Q o T7r){v) >= j[''("'^)'«](^;) (4.6) 

for all C G g. As before, denotes the force field with delay induced by Y . 

We shall see in Section 5 how to construct such a force field with delay 
in the case Q = G. 



5. Dissipation with delay for Euler— Poincare and Lie— Poisson 
equations 

A key step in the reduction of the Euler-Lagrange equations from the 
tangent bundle TG of a Lie group G to its Lie algebra g is to understand how 
to drop the variational principle to the quotient space. The formulation of 
the Euler-Poincare equations and the reduced variational principle is given 
by 

Theorem 5.1. Let G be a Lie group and L : TG ^ IR a left invariant 
Lagrangian. Let I : g ^ ]R be its restriction to the tangent space at the 
identity. For a curve g{t) G G, let^it) = Tg(t)Lg{t)-ig{t). Then the following s 
are equivalent: 

i) g(t) satisfies the Euler-Lagrange equations for L on G. 

ii) The variational principle 

5 f L{g{t),g{t))dt = Q (5.1) 

J a 

holds for variations with fixed endpoints. 
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iii) The Euler-Poincare equations hold, 





(5.2) 



iv) The variational principle 




S l{^{t))dt ^ (5.3) 

J a 

holds on g, using variations of the form — rj + r}] where rj vanishes at 
the endpoints. 

In coordinates, the Euler-Poincare equations are 

= CLq^.C (5.4) 

where C^^ are the structure constants of g relative to a given basis ^" are the 
components of ^ relative to this basis. 

Since the Euler-Lagrange and the Hamilton equations on TQ and T*Q 
are equivalent if the fibre derivative of L is a diffeomorphism from TQ to 
T*Q, it follows that the Lie-Poisson and the Euler-Poincare equations are 
also equivalent under similar hypotheses. To see this directly, we make the 
following Legendre transformation from g to g*: 

dl 

dh 

and assume that £ ^ u is a diffeomorphism. Note that — = E and so it 

is now clear that the Euler-Poincare equations are equivalent to the Lie- 
Poisson equations on g*, namely 

-p = ad*dhii. (5.6) 
dt 

Now we are ready to synthesize our discussions on forces with delay and on 
the Euler-Poincare equations and to transfer this forcing with delay to the 
Lie-Poisson equations by means of the Legendre transform. We begin with 
a formulation of the Lagrange-d'Alembert principle with delay. 
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Theorem 5.2. Let G be a Lie group, L : TG ^ R, a left invariant 
Lagrangian, and F : TG x TG — > T*G a force field with delay equivariant 
relative to the canonical left actions of G on TG x TG and T*G respectively. 
Let / : g — > ]R and / : g x g ^ g* be the restriction of L and F to T^G — g. 
For a curve g{t) G G, let ^{t) — Tg^t^Lg^f^-ig{t). Then the followings are 
equivalent: 

i) g{t) satisfies the Euler-Lagrange equations with forcing with delay for 
L on G. 

ii) The integral Lagrange-d'Alembert principle with delay 

6 j\{g{t),g{t))dt = j'^ ({g{t), g{t)): (mM))) ■ 5g{t)dt (5.7) 

holds for all variations Sg{t) with fixed endpoints. 

iii) The Euler-Poincare equations with forcing with delay are valid 

iv) The variational principle 

5 j'mmt = J J {mm) ■ m)dt (5.9) 

holds on g using variations of the form — f] + [Cv] where r) vanishes at 
the endpoints. 

Proof. We have aheady seen that i) and ii) are equivalent for any configu- 
ration manifold Q in Section 2. Next we prove that ii) and iv) are equivalent. 
First, note that / : g ^ IR and / : g x g — g* determine uniquely a function 
L : TG IR, and a function F : TG x TG — > TG by left translation of 
the argument and conversely. Thus, the equivalence of ii) and iv) comes 
down to proving that all variations Sg{t) e TG of g{t) with fixed endpoints 
induce and are induced by variations S^{t) of ^{t) of the form 5^ — rj + [^, 77] 
where r]{t) vanishes at the endpoints. But this is precisely the matter of the 
proposition 5.1 in [2]. 

The Euler-Poincare equations with forcing with delay have the following 
expression in local coordinates: 
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where Cj„ are the structure constants of the Lie algebra g. 

The condition that the integral curves of the dissipative vector field with 
delay preserve the inverse images of coadjoint orbits by momentum map and 
hence the integral curves of (5.8) preserve the coadjoint orbits of g* is given 
by (4.4), Section 4. Since ^gIq) = T^Rgi^) and J{vg) = T*RgWL{vg), we get 

< F, O TT > {Vg, Vg) =< F (Vg, V g) , T,Rg{^) >= T.^i^.F (v g , V ~g) ' ^ 

and 

JH^^^-^)A{V^) = T:RgWLiVg).[v[Vg,Vg),^] = 

= (^d;^v,^,)°T:R,oWL){vg)-i. 

Since F and WL are equi variant, 

T:RgF{Vg,V~g) = Ad^^ F {Tg L g-,V g , T~g L g-^V ~g) 

and 

° ^^^^ ° = ° ^d;-. O (TgLg-.Vg). 

Because Adg-i o ad^^^^ ~ ^ = ac^^^g-i^^^^,^.) o Ad^-i we get 

jH^^-^U](vg) = o ad*^^^_^^^^^,^^ o FL)(T,L,-..,) 

and the identity (4.4), Section 4 thus becomes 

F (TgLg-iVg,T~L^-,Vg) = {ad*^^ ^^v , ii^) ° (^^^s-^^s) • 

Letting ^ = TgLg-iVg, f = T~L^-iVg, this becomes 

(9/ 

The left hand side is independent of g and thus the right hand side must 
be also gf-independent. Thus taking g = g = e the criterion (4.4), Section 4 
becomes: for every ^, ^ G g there is some rj U,n eg such that 




(5.12) 
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In other words, the force field with delay / (and hence F) is completely 
determined by an arbitrary map ?] : g x g — > g via formula (5.11) and we 
conclude the following 

Corollary 5.3. The solutions of the Euler-Poincare equations with for- 
cing with delay (5.8) perserve the coadjoint orbits of g* provided the force 
field with delay f is given by (5.12) for some smooth map : g x g ^ g. 

Transforming the Euler-Poincare equations with forcing with delay by 
means of the Legendre transformation (5.5), the equations (5.8) with the 
force field with delay (5.12) become 

-4- — ad*ai a — —ad*/ ~\ u (5.13) 

where : g* x g* — > g. The requirement on the map rj is that the right hand 
side of (5.13) be a gradient relative to a certain metric on the orbit. 

To generalize the metric defined by Killing form [3] to coadjoint orbits 
of the dual g* of a general Lie algebra g, we introduce a symmetric positive 
definite bilinear form. 

Let a symmetric positive definite bilinear form F : g* x g* — > IR and 
denote by F : g* — > g the induced map given by r(a,/3) =< /3, Fo; > for 
all a,P e g*, where <, >: g* x g — > IR is the pairing between g* and g. 
Symmetry of F is equivalent to symmetry of F, i.e. F* = F. We introduce 
the following new inner product on g: 

<e,r7>r-i = <r-Ve> (5.14) 

for all ^,ri E g and call it the F~^-mner product. 

Let g/i = ^ g I ad^iJ, = 0} denote the coadjoint isotropy subalgebra 
of /X e g and denote by g'^ the orthogonal complement of g^^^ relative to 
the F~^-inner product. For an element ^ G g we denote by and ^'^ the 
components of C, in the orthogonal direct sum decomposition g = g^ © g^. 
For all g e G, ^ e g and Ad* : g* g*, ao?| : g* g*, Adg : g ^ g we have 

(5.15) 



TiAdln) = AdgVpL, 

adliAdlii) = Adl{ad\,^^{Ad*.,ij)), // e g*. 



Let (5^^^ be the coadjoint orbit through /xq G g* and /j, fi E ^un- There 



exists g E G such that n = Ad* ifl and g^ = Adg-ig~, g^ = Adg-ig^ 
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Let C be a positive Casimir function on g* and let uj be the coadjoint 
orbit symplectic structure defined by 



(5.16) 



for all /X e (5;,^, ^,?7 e g and ?7*(/^) ^^T(T^^. li /i, /i e then 

ad*^li e Tij_&^^, adyx e T~(y^^. Since = Ad*g-i]l it results 

adln = Adl-i (ad\^^^]i) . (5.17) 

We define the (C, V''^)-normal metric on (5^^ respect to /i and Jl by 

< adyi.adyi >n {/J') = C{jl) < T'^r]'^,^'^ > - 



coin) ((r/x)., (r/2),) • ujiTr]) ((r/x),, (r/i),) 



(5.18) 



for all ^, 7/ G g. 

Wc shall regard C and F as fixed in the following discussion and just refer 
to this metric as the normal metric. U n = p, then the normal metric is given 
in [2]. 

Let k : g* — i> IR be a smooth function. We shall compute the gradient 
vector of k\0- with respect to the normal metric. For this purpose we denote 

dk 

by — G g the derivative oi k at /i and by grad/c(//) the gradient of k^o- . 

Since grad/c(//) G T^&^g we can write gradA;(/i) = ad*/! for some rj E g. 
Since and rj^ are orthogonal in the F~^-inner product, we get 



< ad%kii,^ > [jj) = (/X, 



' dfj, 



(n) = /adl^,^\ (/x) 



grad/c(/x),a(i|/x >iv (a*) =< ad*^n,adln >^ (/x) = C(/x) < F ^i'^.r]^ > - 
1 

1 



c{^x) 



< FC, [F/i, F/2] > . < Fry, [F/x, F/I] >= C(/2) < V-\e + ^^), > - 
< F-^[F//,F/i],e > • < F-i[F//,F/i],77 >= C(/I) < F-V,^ > - 



< F-i[F//,F/i],77 > • < F-i[F/x,F/I],e > 
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for any ^ G g. Therefore 

C{]X)V-'7f - -^V-\V^^,VJ^Y < r-^[r/.,r/i],r7 >= -ad> (5.19) 



Because rj'^ — AdgT)'^ it results 



C{Ad*ij)Adgr]^ - ( r ( — ) , [r//, VJl] ) = -r ( ) (5.20) 



and 



or 

= - cm' ("'S'^) ^ co;^ (I) ■ 1^"' ^« > ■ l^"- 

Thus 

gradA;(/x) = ao?*^/^ = - F adafe/I + 

and the equation of the gradient vector field in e relative to the normal 
metric on Cy.,^ is 

^ :</.A/^ + 77r^777-v/rf|^V[r/.,r/i]V[r/.,r/ir (5.23) 



dt C{il) r(g)^^C7(/x)C(/i)\ V^/.^ 



Therefore in (5.23) we put ri{fi, /i) — —r]^ and the Lie-Poisson equations 
with delay forcing (5.6) become 

If g is a compact algebra with the bi-invariant inner product < ■, ■ > on 
g and g is also semisimple, then we could let < •, • >= —k{-,-) where —k{-,-) 
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is the Killing form. In these conditions the inner product identifies g with 
its dual g*, coadjoint orbits with adjoint orbits so that ad^/i = [fJ^,^] and 
dk 

— — Vk{iJ,), where Vk{fi) is the gradient of /c on g at ^ relative to the 

bi-invariant inner product < ■, ■ >. The formula for the gradient vector field 
on the adjoint orbit (5^^ relative to /i and fl becomes 

(5.25) 

where F : g — ^ g defines the symmetric positive definite bilinear form 
{C:V) ^< r^,?7 >. Thus in this case the Lie-Poisson equations with de- 
lay forcing become 

^ = -[Vh{n),n] + ^[/^, r \fi, vkim - 

^ (r(v^(//)), [r/., m) ■ [//, [r//, rjin (5.26) 



c(/i)C(/i) 

Taking C(/i) = 1 and F to be the identity, the dissipative term with delay in 
(5.26) is the Brockett double bracket. The condition that the delay forcing 

term be dissipative is — < and this imposes some conditions on the choice 

at 

of the function /c : g* — > R. 



6. Free rigid body with delay 

Let G = 5*0(3) and the usual identification (so(3), [■, ■] = (IR^, x) implies 
so(3)* = IR'^ via the natural pairing given by the Euclidean inner product. 
Consider (T[Mo) the coadjoint orbit through Mq G 1R^. The infinitesimal 
generator of the coadjoint action is given by Cso{3)* (M) — ^ x M ior M e 
&{Mo) and ^ e so(3). For M,M e ©^(-Mg) let the tangent vectors Cir3(M) = 
^ X M e Tm&{Mq) and ^^^{M) = ^ x M e T^O'{Mq)] the coadjoint orbit 
symplectic structure becomes 

^ (eK3(M),eH3(M)) = M • (m X e) . (6.1) 
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The normal metric on &{Mo) with respect to M and M, given by (5.18), is 
(^M X C,M xr))^{M) = ^{m X (m X i) ,M X (m X ri)) - 

-^^[M-(Mxi))-(M-{Mxri)) = 

^{mx^,Mx n)^ (m) - i (m . (m X e)) • (m . (m X n)) (6.2) 

with c = C(M) = C(M). _ _ _ 

The normal metric at M on two tangent vectors M x ^, M x 77 to the 
sphere of radius c is given by 

(^M x^,M xri)^(^M) = ^(^M X (m X ,M X (m X f])) (6.3) 

where the inner product of the right hand side is the standard inner product 

in IRl _ _ 

The normal metric at M with M 7^ M, on two tangent vectors M x ^, 
M X ?7 to the sphere of radius c at M is given by 

(m X ^, M X T])^ (M) = -J (m X (m X ^) , M X (m X 77)) (6.4) 

where the inner product of the right hand side is the standard inner product 

in B?. From (6.3) and (6.4) it results (6.2). 

1 

From (5.26) where k{M) = h{M) = -\\Mf it results the Lie-Poisson 
equation for the rigid body with delay: 

M{t) = M X Q + X (m X Q) . (6.5) 

We shall discuss the stability of the equilibrium states for a free rigid 
body with delay. 

Let the free rigid body with delay given by the equation 

M ^ M xn + aM X (^M xfl^ , (6.6) 

where M = m= {hx{t), hvit), l3z{t)f , O = {x{t),y{t), z{t)f , M = m = 
{hx{t - r), hyit - t), hz{t - t)Y, n = {x{t - T),y{t - r), z{t - T)f, h > 
0, /2 > 0, /a > 0, a a constant and r > 0. 
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It is not hard to see that the equihbrium states of our system are Oi = 
(m//i,0,0)^, ^2 = (0,m//2,0)^, = {0,0,m/h), m G W. 

Proposition 6.1. The equilibrium state fti has the following behavior: 
(i) The corresponding linear system is given by 



5n = ASn + aG5VL 



(6.7) 



where 





( ° 
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( ° 
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Ills 


J 



(6.8) 



A 



ii) The characteristic equation is 
ani^ fh — h , h — h 



A 



h 



+ 



Ae 



-rA 



+ 



2 4 

am 



(/2-/i)(/3-/i)e- 



-2rA 



(/l-/2)(J3-/l) 2 

— m 



= 0; 



(6.9) 



Hhh 

(iii) On t/ie tangent space at VLi to the sphere of radius rn^ the linear 
operator given by the linearized vector field has the characteristic equation 

.2 ^i^-i^ i^-i^- 



am 



+ 



Ae 



-tA 



+ 



2 4 

am 
Ilhh 



:/2-/i)(/3-/i)e-'^^- 



(/l-/2)(/3-/l) 2 



0; 



(iv) // /i > I2, h > h for < r < Tc, where 

h [h{h-h) + h{h-h) 



2>\a\m^{h-h){h-h) 



(6.10) 



(6.11) 



then the equilibrium state Q^i is asymptotically stable. 

Proof, (i), (ii), (iii) result from the definitions of the liniarized and the 
characteristic equation by calculus. For (iv) consider Ii > I2, h > I3. If 
r = then the characteristic equation (6.9) has eigenvalues with the real 
parts negative and fli is asymptotically stable. Following [9] and [12] it 
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results that for < r < Tc the equihbrium state fli remains asymptotically 
stable. 

In the following we study the existence of Hopf bifurcations for the free 
rigid body with delay (6.6) by choosing the delay r as a bifurcation para- 
meter. First we would like to know when the equation (6.10) has purely 
imaginary roots ±ia;o (ti^o > 0) at r = Tq. Note that A = icuo is a root of (6.1) 
if 

u^-c- auo sin uoTq - b cos 2a;oTo = 0, 

(d.Iz) 

owq cos ujqTo — b sin 2ujoTo — 



2 2 4 

am , ^ , , , a m 



with a, 6, c given by 

t2t t L-'Sl-'i -zj I -ayj, j2t j 

-'l-'2-'3 Jl-'2-'3 

2 

m^ 
Ilhh 



c^-^{h-h){h-h). (6.13) 



We deduce the following 

Proposition 6.2. (i) If\m\ < - — r then A = iujq is a simple root o/(6.10) 

\a\ 



and 



a + ^o?-A{b-c) TT 
"^o = ^ , To = — , To > Tc. (6.14) 

(ii) // \m\ > - — - then A = icoo is a simple root of (6.10) and 

\a\ 



-a + ^a'^-^b-c) 37r ,^ ^ , 

(^0 = ^ , To = — , To > Tc. (6.15) 

We proceed to calculate Re ( — ] at r = tq. By differentiating the 

\dTj 

equation (6.10) implicitly with respect to r, we obtain: 

dX _ aX^e-^^ + 26Ae-^^^ 

~ 2A + a(l - rA)e-^^ - 26x6-^^^ ' ' 

It is then evaluated at A = iujq and t — tq given by (6.14), (6.15), yielding 

Re f = a;o(^o + a)(a -26) 
\ar J x=iu,o Tolauo — 2by + [uo + ay 

T=TQ 
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From the standard Hopf bifurcation theory we have the following result. 
Proposition 6.3. If Ii > I2, Ii > I3 and ujq, tq are given by (6.14), (6.15) 

with (6.13), then Re ( — ] 7^ and a Hopf-type bifurcation occurs at 



dr I x=iu. 



Jli when r passes through tq. 

In the following we obtain some conditions which guarantee that the free 
rigid body with delay undergoes a Hopf bifurcation at r = tq. The method 
we use is based on the normal form theory and the center manifold theorem 
introduced in [7]. 

With the translation V — Q—Qi, U — M—Mi the equation (6.6) becomes 
U = lAV + alGV + F (u, V, U, v) , (6.18) 
where A, G are given by (6.8), / = diag (/i, h) and 

F (U,V,U,V) =U xV + a[u X (U xVli) + U X (Mi xV) + 

+ Mix(Ux V)] +aU X (U xV) . (6.19) 
Prom (6.18) it results 

V = AV + aGV + N (v, v) , (6.20) 

where N (v, v) = r^F (/- V, V, I'^V, v) . 

For (f) e C^([— To, 0], IR^) we define an operator A by 

r d(f) 



dO^ (6.21) 



{ A(t>{Q)+aG(t>{-To), e^O 
and for ijj e C^([0, Tq], H^*) we define the adjoint operator A* of A by 

"dJ' '^^^'^'^ (6.22) 

i^{0)A + aij{To)G, s = To; 

A and A* are adjoint operators with respect to the bilinear form 

< V, >= mm - « r r ^(^ - o)Gmd^d9, (6.23) 

J -TO J^=0 
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e C\[-ro,0],¥i'), i' e C\[0,TolM^*)- 

Let now Ai = iuiQ, A2 = Ai = —iujQ be eigenvalues of A, where o^o is given 
by (6.14). They are also eigenvalues of A*. We can easily obtain that 

m = (0, V2, vsfe^^', e e [-To, 0], (6.24) 

where V2 — {I3 — /i)"^, = Xihh — {h — h)m'^ae~^^'^° , is an eigenvector of 
A corresponding to Ai and 

jP{s) = (0, W2, ws)e^'\ s e [0, To], (6.25) 

where W2 = hih—h)^, = (Ai/1/2 — (/2— /i)m^ci;e^^^")/3, is an eigenvector 
of A* corresponding to Ai. 

Prom (6.23), (6.24), (6.25) it results 

On = <W:(P>^'^2W2 + V3W3--^— rri^^ + 

+ Ase^^ro _ _ I^)V2W2 + hih - h)v3W3i 

2 

ai2 = <V,^>='y2«^2 + ^^3«^3--^^^ -2^(2-6-^2^°- (6 26) 
- e^^^^)[h{l2 - h)v2W2 + hih - hyo^w^l 

021 = < V^, >= < V', > = ai2, 022 =< i>-, > = 



= <V',0> = aii. 

Let d — aiiOii — 012012 and bu — 612 = — The vector field given 

da 

by tjj(s) = bui'is) + 6i2'?/'(-s), s e [0,ro] is an eigenvector of A* satisfying the 
relations: 

<tp,e>=i, <i>,e >=<i),d <i>,e>=i. (6.27) 

Then we have 

i^{s) = (0, W2, W3) e^^^ s e [0, To], (6.28) 

where 

W2 = {bu + 612)^2, W3 = 611W3 + 612W3. (6.29) 

We shall discuss the existence of a local center manifold around the equi- 
librium point of the equation (6.20). 
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Let B = CH[-T-o,0],]R^) and A the vectorial space, span of the eigen- 
vectors (f){9), (t){0) corresponding to Ai = iuq, A2 = —iujQ. For a given 
neighborhood V of e B, a local center manifold W^^^{Qi) — W'^{Q, V) of the 
equilibrium point 0(0, 0, 0) of (6.20) is a C^-submanifold that is a graph over 
V n A in i3, tangent to A at O and locally invariant under the flow defined 
by the equation (6.20). In other words 

W^^iQi) ^{ip^ eB\ip'' ^u(t) + u4> + w{ip), e V n a} , (6.30) 



where w : A — > B is a C-'^-mapping with (/?(0) = 0, ^^^^^(O) = and {^l), wj — 

0. Moreover, every orbit that begins on W,oc(0) remains in this set as long 

as it stays in V. 

The basic result on the existence of the local center manifold for the delay 
differential equations is given in [6] . From the definition of the local center 
manifold it results that 

W^ioc(O) n Vi = ^Lp" ^ B \ if" ^ u(f) + u(t) + w{u, u), U^Ui + iU2, 

{UUU2) eVi CH^} , (6.31) 



where w : [— t, 0] x (D IR is given by w{6,u,u) = w [u(f){6) + u 
Following [6], ip'^ G W^ioc(O) H Vi and the function w is the solution of the 
partial derivate system 

^ {9Mt)Mt))+9{u{t),u{t)) m +^WhW})-m = ^ {9Mt)Mt)) 

(6.32) 

with 

dw 



(0, u{t),u{t)) + g {u{t),u{t)) m + 9 {u{t),u{t)) ■ 0(0) 



= Aw (0, u{t),u{t)) + aGw (-tq, u{t),u{t)) + N {u{t)(t){0) + u{t)(f){0) + 

+ w (0, u{t),u{t)) , u{t)(t){-TQ) + u{t)4>{-TQ) + w (-To, u{t),u{t))) , (6.33) 
where u{t) is a solution of the ordinary differential equation 

u{t) = Xiu{t) + g {u{t),u{t)) (6.34) 
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and 

g {u, u) = i^{0)N (u(j){0) + 110(0) + w (0, u, u) , u(f){-To) + u'${-to) + 

+ w {-To,u,u)) . (6.35) 
For ip'^ G W^ioc(O) n Vi the solution of (6.20) is given by 

v{t){e)^u{t)(f>{e) + u{t)^{e) + w{e,u{t),u{t)), oe[-To,o]. (6.36) 

Because the equation (6.20) has a Casimir function (conservation laws) = 
m. for the solution given by (6.36) we have IV{t) ■ IV{t) — 0. 
Consider the function w given by 

w {9, u, u) = -W2o{9)u^ + wii{9)uu + -wq2{9)u^ ■, (6.37) 

with wo2{9) = W2o(0), wii(^) = wii{9), 9 G [-ro,0]. From (6.20) it results 
that the components of ^V^, for V = {x^, x'^, x^)^, V — {x^, x'^, x^)^ are 
the foUowings: 



h- h 2 3 , 

-X X -\- am 



hih - -^2)^2-2 _ /3(-^3 - -^l) ^3~3 



+ 



h{h — h) 1 2-2 -^3(-^3 - A)~i 3-3 

Ju Ju Ju ^ Ju Ju Ju 

h 



h- h 13, 

— X X + am 

2 



+a 



h 

N' iv, V) = ^-^xW + 
-\-a 



XXX — 



am 



XXX — 



h 




hih - 


h) 


h 




h{h - 


h) 


h 




h{h - 


h) 


h 




Hh - 


h) 



XXX 



(6.38) 



XXX 



From (6.38) with V{t), V{t) = V{t){-TQ), given by (6.36) and (t){9), (f){9) 
given by (6.24), it results 

N (V{t), V{t)) = ^F2o + Fiiu{t)u{t) + ^Fo2u{tf + ^F2iu{tfu{t) , (6.39) 



29 



with F20 - (-^20) -^2^0' -^2^o)^> -^11 — (-^11) -^11) -^n)^) -^02 - (-^02) -^0^2) -^0^2)^) 
F2i = {Fi„Fl,Fi,f, where 



-'1 J-1 
-^2^0 ~ = 0, 



-^11 — 



h-h 
h 



am 



V2 (^3 + ^"3) + h{h - h)V' 



Fh - Ff, = 0, 



c,i 2(72 - -^3) _ , c^m 
F02 = ; ^^2^3 + 



h 

3A1TO _j_ gA2T0^ 



h 

p2 _ p3 _ Q 



h 



l2{h-h)vl-h{h-h)vl 



^Airo 



i^2\ = 



+ 



/2-/3 
h 



+ 



2aml2{h-l2) 
h 

2amh{h-h) 
h 



(6.40) 



Fix = —J ^^3W^o(0) -J W2^o(0)e^^^°, 

-'2 -'2 

-'3 -'3 

From (6.35) with '^(0) given by (6.28) andiV(v^(t),1/(t)) given by (6.39) 
we obtain 

g{u{i),u{i)) = ]^g2xn{ifu{i), (6.41) 

where 

92x^W2Fl^^^3Flx. (6.42) 

Taking into account of (6.32) it results that W2o{0), Wn{0) verify the 
differential equations 



w2o{e) = 2XxW2{e), wxi{e) = o, ee [-ro,o]. 



(6.43) 
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From (6.43), (6.33) and IV {t) ■ IV {t) = it results 

W2oie) = E,e^^'', wu{0) = 0, ee [-To, 0], (6.44) 
where Ei is the solution of the linear system of equations 

(a + ae-^'^'G - 2XiE) Ei = -F20. (6.45) 
From (6.45) and (6.40) we deduce 

^oW = ^^2V^', W2oW=^2oW = 0, ^e[-To,0], (6.46) 

obtaining the following 

Proposition 6.4. The solution of the equation (6.6) near upon the sta- 

— ,0, Oj is 

x{t)^- + Re{w'2oiOyit)), 
-'1 

y{t) ^2v2Re {u{t)\ (6-47) 
z{t) = 2Re{v3u{t)), 

where u{t) is a solution of the equation 

u{t) = ^g2iuitfu{t) (6.48) 

and W2q{0) = 777— '^2,V3 are given by (6.24), g2i is given by (6.42). 

Based on the above analysis and calculi, we can see that g2i, wl(^{0), V2, V3 
are determined by the parameters and the delay of (6.6). Thus we can 
explicitly compute the following quantities: 

r(o)-\ u - ^^(^^(Q)) 



X=iu}Q 
r=TQ 



/m(Ci(0))+/X2/m(f),_^ 
T2 = P2 = 2Re(Ci(0)). 
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In summary, this leads to the following result. 

Proposition 6.5. In the formulas (6.49) /X2 determines the direction of 
the Hopf bifurcation: if 112 > ( respectively 112 < 0) then the Hopf bifurcation 
is supercritical (respectively subcritical) and the bifurcating periodic solutions 
exist for r > Tq (respectively r < Tq); (32 determines the stability of the 
bifurcation periodic solutions: the solutions are orbitally stable (respectively 
instable) if 132 < Q (respectively P2 > 0^; T2 determines the periods of the 
bifurcating periodic solutions: the periods increase (respectively decrease) if 
T2 > (respectively T2 < 0). For h = 0.8, I2 = 0.5, I3 = 0.4, a = 0.3, 
m = 1.5 and uq, tq given by the formulas (6.14) we obtain ujq = 3.20631, 
To = 0.88154, /X2 = 0.00958, T2 = 0.00057, (32 = -0.00139. The limit cycle is 
supercritical with the period T2. For Ji = 0.8, I2 = 0.5, /a = 0.4, a = 0.3, 
m — 1.8 and cuq, Tq given by the formulas (6.15) we obtain cuq — 0.68547, 
To = 0.88154, /i2 = 0.00344, T2 = 0.00050, P2 = 0.00097. The limit cycle is 
supercritical with the period T2. 

7. Conclusions and comments 

In this paper we have given a general method to construct a dissipative 
mechanism with delay preserving the symplectic leaves of the reduced space 
and dissipating the energy. The most important case is that of the dual of a 
Lie algebra when the dissipative term with delay is shown to have a double 
bracket with delay form. This theory applies to a number of interesting 
examples from ferromagnetics, ideal fluid flow and plasma dynamics in which 
the previous state of the phenomenon is important. In the future we would 
like to analyze the systems given by Landau-Lifschitz equations with delay, 
the Heavy Top with delay etc. 
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